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In todays rotordynamic calculations, the input parameters for a Finite Element Analysis 
(FEA) determine very much the reliability of Eigenvalue and Eigenmode predictions. While 
modelling of an elastic structure by means of beam elements etc. is relatively straightfor- 
ward to perform and the input data for journal bearings are usually known exactly enough, 
the determination of stiffness and damping for labyrinth seals is still \ the subject of many 
investigations. Therefore, the rotordynamic influence of labyrinths is often not included 
in FEA for rotating machinery because of a lack of computer programs to calculate these 
parameters. This circumstance can give rise to severe vibration problems especially for 
high performance turbines or compressors, resulting in remarkable economic losses. 

The forces generated in labyrinths can be described for small motions around the seal 
center with the linearized force-motion relationship of eq. 1: 
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where K, k are the direct and cross-coupled stiffness and D, d are the direct and cross- 
coupled damping. These elements are the so-called dynamic labyrinth coefficients. 

Several years ago, we started with the development of computer codes for the determina- 
tion of rotordynamic seal coefficients. The following sections will introduce our 
different approaches to evaluate the dynamic fluid forces generated by turbulent, com- 
pressible labyrinth seal flow. 
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NOMENCLATURE 
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Direct and cross-coupled stiffness 

Direct and cross-coupled damping 

Rotor displacements 

Rotor velocities 

Labyrinth fluid forces 

General variable 

Diffusion coefficient 

Source term 

Cylinder coordinates 

Transformation coordinate 

Perturbation parameter 

Rotor eccentricity 

Rotor precession frequency 

Shaft rotational frequency 

Local seal clearance 

Nominal clearance 


G Labyrinth leakage 

r Radius from seal center 

a 

to stator 

r. Radius from seal center 

to end of seal fin 
x,y cart, coordinates 

x,y rot. coordinate system 

Wj preswirl 

t shear stress 

u,v,w axial, radial and circum- 

ferential velocities 


Indices 

0 Zeroth order 

1 First order 


a before labyrinth 

b behind labyrinth 


COMPUTATIONAL FLUID DYNAMICS METHODS (CFD) 

The compressible, turbulent, time dependent and threedimensional flow in a labyrinth seal 
can be described by the Navier-Stokes equations in conjunction with a turbulence model. 
Additionally, equations for mass and energy conservation and an equation of state are 
required. These equations can be solved with the aid of a finite difference procedure. To 
determine the desired coefficients, two methods have been developed: 

- 2-dimensional procedure based on a perturbation analysis 

- 3-dimensional theory using a moving frame of reference rotating with the shaft. 
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PERTURBATION ANALYSIS 


To describe the labyrinth seal flow, we use the time-averaged conservation equations for 
momentum, mass and energy and the equation of state for a perfect gas. The correlation 
terms of the turbulent fluctuation quantities are modeled via the k-s turbulence model 
of LAUNDER and SPALDING /l/. All these equations can be arranged in the following 
generalized form: 
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where 0 stands for any of the dependent variables(e.g. velocities, pressure, temperature, 
density), is the diffusion coefficient and S^ the source term. To calculate the dynamic 
coefficients, several assumptions are introduced: 

1 . small rotor motions around the seal center on a circular orbit 

2 . coordinate transformation 

3. perturbation series expansion for the dependent variables 

4 . solutions for the dependent variables corresponding to the temporal and circum- 
ferential variation of the seal clearance function 

Assumption 1 assures that we can use eq. 1 to calculate the dynamic coefficients. 

When supposing that the shaft rotates on a circular orbit, one can imagine that the local 
seal clearance h varies in time for every location <p=const. To avoid the use of time 
dependent calculation grids, we transform the governing equations to another coordinate 
system by introducing a new radial coordinate 7] whereby the eccentrically moving shaft 
is converted to a shaft rotating in the center of the seal (see Fig. 1). 

The third assumption, together with assumption 1, implies that the perturbation series 

<D = e 0 ®^ + e 1 0. + e 2 + ... + e n <E> n = 1 ... cc (3) 

O 1 Z n 

for the dependent variables can be truncated after the linear term. 
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Fig. 1: Coordinate transformation 

Inserting these equations for every dependent variable into the basic conservation relations 
yields two sets of equations. The zeroth order set describes the flow in the centric laby- 
rinth while the first order equations are governing the flow field for small eccentric shaft 
motions. The last assumption allows the elimination of thei temporali and circumferentialideriva 
tives in the first order equations analytically. 

The seal clearance function for a circular shaft precession orbit around the seal center 
can be stated as 

h = h Q + eh 1 = C r - X(t)cos<p - Y(t) sin<p (4) 

Corresponding to this function, we assume that the dependent variables vary in the same 
way concerning the circumferential direction. 

Thereby, the circumferential derivatives can be calculated analytically. After separating 
the equations into sine and cosine terms and rearranging them by introducing complex 
variables, solutions for the first order variables corresponding to the temporal change 
of seal clearance due to a circular shaft orbit are prescribed, sallowing 1 the eliminat ! ion:of the 
derivatives with respect to time. Finally, we obtain two sets of differential equations which 
have to be solved subquently in order to determine the desired coefficients. 


THREEDIMENSIONAL METHOD 

In contrast to the perturbation analysis, where we have made many assumptions to reduce 
the computational effort, the threedimensional procedure shows much greater generality 
in avoiding these restrictions. Again, the time-averaged conservation equations are solved 
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in conjunction with the already mentioned k-s model. 

To determine the dynamic coefficients we assume that the shaft moves on a circular orbit 
with precession frequency Q around the seal center. Since this would normally result in 
a time dependent problem we introduce a rotating coordinate system which is fixed at 
the shaft center (Fig. 2). In this moving frame of reference, the flow is stationary. Due 
to the rotating coordinate system, centrifugal and coriolis forces occur in the equations 
for radial and circumferential momentum. 


y 4 



Fig. 2: Rotating coordinate system 


SOLUTION PROCEDURE FOR THE CFD-METHODS 

Because of the complexity of the derived equations, no analytic solution can be obtained. 
Therefore, there is a need for a numerical algorithm to solve the equations with the aid 
of a computer. We made our decision for the Finite Difference Method (FDM), which is 
well established in fluid dynamics. 

The FDM procedure starts with the discretization of the calculation domain, which must 
be performed in a different manner for the two methods: While we can use a twodimen- 
sional Finite Difference grid for the perturbation analysis (Fig. 3), the threedimensional 
approach requires additional grids in the circumferential direction (Fig. 4). 


207 




Fig. 3: Grid for perturbation method 



Fig. 4: 3D Finite Difference grid 



In the next step, the generalized equation (2) is integrated over a control volume linking 
every node point with his neighbours in space. The resulting algebraic equations are then 
solved for the zeroth and first order equations (perturbation analysis) or the 3D equations 
respectively. 

Finally, the dynamic forces are calculated from the pressure distribution on the rotor sur- 
face for two precession frequencies 0=0 and Q=g), resulting in two sets of forces. Then 
the dynamic coefficients are determined from eq. 1, where the prescribed circular shaft 
orbit is introduced: 


- Fj = r 0 ( K + 0 d ) 

- F 2 = r 0 ( - k + O D ) 


(5) 
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For more detailed information concerning the perturbation analysis or the 3D algorithm, 
readers are referred to WEISER and NORDMANN /2,3/. 


SIMPLIFIED METHODS 

While the CFD algorithms require a lot of computer storage and cpu time, the simpli- 
fied methods described in this chapter show a reduced effort to determine the rotordyna- 
mic labyrinth coefficients. 


THREE VOLUME BULK FLOW MODEL 

Following the conventional approach (see for example IWATSUBO /4/ or CHILDS /5/), 
a three volume bulk flow model has been developed based on the conservation equations 
for mass, momentum and energy working with velocities, pressure, density and tempe- 
rature, which are averaged over the control volume height. 

Flow visualization experiments show, that for look-through labyrinths, the flow field can 
be divided in two characteristic regions: a vortex flow in the seal chamber and a jet flow 
region beneath the seal strip and the groove. Therefore we use three control volumes (CV) 
to describe the flow situation in the labyrinth (Fig. 5). 

To account for the flow turbulence, wall shear stresses and a fluid shear stress in the 
contact region of CV II and CV III are introduced. Following SCHARRER /6/, a vortex 
velocity in the chamber is assumed. The constants appearing in the shear stress formula- 
tions are determined with the aid of the CFD programs for centric shaft position. 

Again, a perturbation analysis in connection with the implementation of solutions for the 
temporal and circumferential variation of the dependent variables is performed, resulting 
in the already known two sets of differential equations, which now depend only on the 
axial coordinate z. After their numerical solution, the forces and the dynamic seal coeffi- 
cients are calculated by a pressure integration over the shaft surface. 

More information concerning the bulk flow model can be found in NORDMANN and WEISER 
77,8/ . 
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Coordinate System 
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Fig. 5: Three volume bulk flow model 


CORRELATION EQUATIONS FOR THE DYNAMIC COEFFICIENTS 

Already in the design stage of a turbine or a compressor, the estimation of labyrinth 
coefficients is important for inclusion of these effects into rotordynamic calculations. 
Therefore, we developed approximation formulas for the desired stiffnesses and dampings 
to provide a tool for the design engineer, which can give some information about the 
rotordynamic seal influence to be expected. These correlations are based on extensive FD 
calculations for a model labyrinth, where performance data and geometry were varied in 
a wide range of practical interest. 

First, the influence parameters which will be included in the correlation to perform must 
be defined (shown in Fig. 6). 
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Fig. 6: Influence parameters for the correlation of dynamic coefficients 


For the approximation equations, several definitions are used: 
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DESCRIPTION OF DIRECT STIFFNESS 


The approximation for the direct stiffness consists of three parts: 

- stiffness for axial flow situation ( K 0 ) 

- stiffness induced by shaft rotation ( K^) 

- stiffness induced by preswirl velocity ( K ) 

Wj 

From FD calculations, we found that the influence of preswirl changes sign for labyrinths 
with more than 5 strips and that K becomes negative for multi-chamber seals. This is 
included in the following formulas: 


with 
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CROSS-COUPLED STIFFNESS 

The cross-coupled stiffness is directly related to the leakage loss of the labyrinth. Also, 
the influence of geometry, shaft rotation and preswirl are included: 
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( 8 ) 
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The average circumferential velocity w m is determined using the 1/7-power law. 


DIRECT DAMPING 


To describe the direct damping D, some of the previously defined functions are used. The 
influence of shaft rotation, preswirl and geometry are found to be smaller than for the 
stiffnesses: 
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The correlation coefficients were determined from the FD calculation results for the data 
variations of the chosen model labyrinth: 


C* = 3.141 10" 2 C* = 16.0 Kg/( m s ) C* 

q* _f 4.2 10 -3 ( strips on rotor ) * 

4 ^ 3.5 10 -3 ( strips on stator ) 5 

C* = 0.18 Ns 2 / m 2 C* = 1.60 10" 2 


= 2.80 10" 2 Kg 

1.50 m" 1 C* = 1.571 10~ 3 

N s 2 / m 
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COMPARISON OF THE DIFFERENT APPROACHES WITH MEASUREMENTS 

1. Example: A three chamber look-through labyrinth 

To compare the two FDM, the bulk flow model and the correlation equations to measure- 
ments, we have chosen a look-through labyrinth which was experimentally investigated 
by BENCKERT 79/ (seal data given in /2/). 
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Fig. 7 : Direct Stiffness 



□ 
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Fig. 8 : Cross-coupled Stiffness 


2. Example: A multi-chamber look-through labyrinth seal 

SCHMIED /10/ published experimental results obtained by CHILDS for a 12 strip stator 
labyrinth (seal data given in /10/) . The comparisons in Fig. 9-11 show the results of 
perturbation analysis, bulk flow model and correlation for the stiffnesses and the direct 
damping. 
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Fig. 9 : Direct Stiffness 





APPLICATION OF SEAL DYNAMICS CALCULATION TO ROTORDYNAMIC DESIGN 

AND VIBRATION INVESTIGATIONS 


In 1988, SCHMIED /10/ published rotordynamic calculations for a high performance injec- 
tion compressor, where he included the most important fluid interaction influences coming 
from journal bearings, oil ring seals and labyrinths. Fig. 12 shows a sketch of the compressor 
rotor discretization and gives the basic data. The labyrinth data are given in /10/. SCHMIED 
uses the program MAD YN to perform the eigenvalue calculations. In this chapter, a compa- 
rison is made to SCHMIEDs results for the first forward bending mode with straight 
through shroud and hub labyrinths, the comb-grooved balance piston and ideally floating 
oil ring seals. 
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Basic compressor data: 

Suction pressure: 198 bar 

Discharge pressure: 700 bar 

Mol. weight of gas: 20.05 

Fig. 12: Investigated rotor 
COMPARISON TO SCHMIED'S RESULTS FOR THE SIXTH STAGE SHROUD LABYRINTH 


Mass flow: 
Rotor mass: 
Bearing dist.: 


47.5 kg/s 
26S kg 
1565 mm 


Op. speed: 13400 rpm 


SCHMIED uses the labyrinth seal model of WYSSMANN /II / to determine the rotordynamic 
coefficients. For the sixth stage shroud labyrinth, his calculations show negative direct 
stiffness values. As can be recognized in Fig. 13 to 15, the results of our investigations 
with the perturbation analysis, the bulk flow model and the correlation equations yield 
positive direct stiffness as expected. The agreement of cross-coupled stiffness and direct 
damping to the values obtained with WYSSMANNs theory is good. Especially, the different 
procedures presented in this paper are in relatively good agreement with each other for 
all rotordynamic parameters. This statement holds also true for all other labyrinth seals 
of the investigated compressor. 
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Fig. 13: Direct Stiffness 
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Fig. 14: Cross-coupled Stiffness 
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Fig. 15: Direct Damping 
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ROTORDYNAMIC ANALYSIS OF EIGENPROBLEM AND STABILITY 


SCHMIED investigated three cases of rotordynamic interest: 

- "ideal": rotor running in journal bearings, no labyrinth influence 

- "with swirl brakes": labyrinth influence included for low preswirl conditions 

- "without swirl brakes": labyrinth influence included for high preswirl conditions 

For these cases, we have calculated the eigenfrequencies and modal dampings of the rotor 
system using the Finite-Element program of DIEWALD /12/ . In Fig. 16 and 17, comparisons 
are presented for the first forward bending mode. For the "ideal" system, the calculations 
show nearly the same results. Because of the labyrinth coefficients obtained with the FD 
perturbation analysis for the hub and shroud labyrinths and the bulk flow model for the 
balance piston, the first eigenfrequency is higher than in SCHMIEDs calculations. This 
is mainly due to the positive direct stiffnesses for the impeller labyrinths. Also the modal 
damping shows great differences. 
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Fig. 16: Comparison of first eigenfrequency 
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Fig. 17: Modal damping of first forward bending mode 


CONCLUSION 

Four calculation methods have been described allowing calculation of the dynamic labyrinth 
coefficients. Comparisons to experiments show good agreement for all methods and all 
investigated cases. 

The practical example of rotordynamic calculations for a high-performance compressor 
has shown that the determination of labyrinth coefficients must be performed as exactly 
as possible in order to predict the eigenfrequencies and modal dampings correctly. Also, 
it has been proved that the methods presented in this paper are superior for this task 
to the existing theories. 
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